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Recent studies have shown that concepts of effective field theory such as nat- 
uralness can be profitably applied to relativistic mean-field models of nuclei. 
Here the analysis by Friar, Madland, and Lynn of naturalness in a relativistic 
point-coupling model is extended. Fits to experimental nuclear data sup- 
port naive dimensional analysis as a useful principle and imply a mean-field 
[ expansion analogous to that found for mean-field meson models. 
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I. INTRODUCTION 

Quantum Chromodynamics (QCD) is believed to be the underlying theory of hadrons 
and their interactions and therefore of nuclei as well. However, direct solutions of QCD for 
nuclei (e.g., lattice calculations) are not presently feasible. On the other hand, the modern 
perspective of effective field theory (EFT) [|l]-0 provides a framework for constructing theo- 
ries based on observed hadronic degrees of freedom that can faithfully reproduce low-energy 
QCD. An example of a successful EFT is chiral perturbation theory (ChPT) which 
systematically describes low-energy hadronic processes in the vacuum. 

The success of effective chiral lagrangians in the vacuum has motivated the application 
of EFT concepts to models of nuclear properties. A chiral effective lagrangian for nuclei 
was presented in Ref . M , and has lead to new insights into the successes of quantum hadro- 



dynamic (QHD) models JTOjll] of nuclei. In that work, a general model in which nucleons 
interact via mesonic degrees of freedom was constructed and applied to nuclei at the one- 
baryon-loop level (Hartree approximation). 

The effective lagrangian and energy functional were organized according to Georgi's naive 
dimensional analysis (NDA) which provides a controlled mean-field expansion in 



terms of powers and derivatives of the mean fields if the dimensionless coefficients identified 
in the NDA are of order one. The latter property is called "naturalness" in this context. 
Detailed fits to nuclear observables validated this expansion and truncation scheme f9|. This 
result is somewhat surprising since the NDA counting is based on absorbing short- distance 
physics into the coefficients of the effective lagrangian, but a mean-field functional fit to 
nuclear data must also absorb long-distance many-body effects. 
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A useful interpretation of the finite-density mean-field effective theory is in terms of 
density functional theory (DFT) JlJ -17]. In a DFT formulation of the relativistic many-body 
problem, one works with an energy functional of scalar densities and vector four- currents. 
Minimization of the functional gives rise to variational equations that determine the ground- 
state densities. By introducing a complete set of Dirac wave functions, one can recast these 
variational equations as Dirac equations for occupied orbitals; the single-particle hamiltonian 
contains local scalar and vector potentials, not only in the Hartree approximation, but in 
the general case as well. The Hartree approximation only limits the form of the potentials. 
Thus the effects of many-body physics such as short-range correlations are incorporated 
(approximately) by a direct fit to nuclear observables. 

In mean-field meson models, the scalar and vector meson fields play the role of auxiliary 
Kohn-Sham potentials. Alternatively, one can expand the local potentials directly in terms 
of nucleon densities. The corresponding lagrangian replaces meson exchange with point- 
coupling (contact) interactions. In Ref. [0, such a relativistic point-coupling lagrangian 
was introduced and applied to nuclei at the mean-field level. Subsequently, Friar et al. [19| 



examined this model for naturalness and concluded that the parameters were, in fact, mostly 
natural. Here, the analysis of point-coupling models in the context of EFT's is reexamined 
and broadened in a manner consistent with the analysis of the chiral mean-field meson model 
of Ref. 0. 

We find that the assumption of naturalness is justified and that the density expansion 
implied by NDA is applicable. In principle, the parameters found here could be related to 
those in the meson model of Ref. by using the equations of motion in the lagrangian 
of || to systematically eliminate the meson fields. However, the point-coupling parame- 
ters are underdetermined by the fits, which limits the usefulness of the comparisons. New 
optimization procedures described below suggest that a different organization that exploits 
cancellations characteristic of relativistic models may be more productive. 



II. THE POINT-COUPLING LAGRANGIAN 

At present we cannot derive an effective point-coupling lagrangian directly from QCD. In 
the same spirit as ChPT and the effective meson model of Ref. |§, a general point- coupling 
effective lagrangian is therefore constructed consistent with the underlying symmetries of 
QCD (e.g., Lorentz covariance, gauge invariance, and chiral symmetry). In this work, we 
construct a one-loop energy-functional from the lagrangian and determine the couplings 



by fits to nuclear observables. As noted above and discussed in Ref. |pT| , this approach 
approximates a density functional that, if sufficiently general, incorporates many-body effects 
beyond the Hartree approximation. We expect this approximation to be reasonable because 
of the large scalar and vector potentials ("Hartree dominance") P,pT[]. 

This framework will only be useful if we can identify a valid expansion and truncation 
scheme. This requires an organization of terms in the effective lagrangian and a way to 
estimate the couplings. While precise relations between these couplings and the underly- 
ing QCD parameters are unknown, an estimate of the magnitude of the couplings can be 
obtained by applying Georgi's naive dimensional analysis (NDA) |TJ||5|.|2(| . 

The procedure is to extract from each term in the lagrangian the dependence on two 
primary physical scales of the effective theory, the pion decay constant, f n ~ 94MeV and 
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a larger mass scale, A ~ An f w / ' \J~N~f (where Nj is the number of light flavors) flT2yi3| . The 
scale A is associated with the mass scale of physics beyond the Goldstone bosons (pions): 
the non-Goldstone boson masses or the nucleon mass. This mass scale ranges from the scalar 
mass ~ 500 MeV to the baryon mass ~ 1 GeV. When a specific value is needed in this work, 
A will be taken to equal the p-meson mass (770 MeV), roughly in the center of this range. 
To establish the canonical normalization of the strongly interacting fields, an inverse factor 
of f n is included for each field and an overall factor of / 2 A 2 fixes the normalization of the 
lagrangian. The physics of NDA is discussed further in Refs. and [11|. 



We construct the point-coupling effective lagrangian as an expansion in powers of the nu- 
cleon scalar, vector, isovector-vector, tensor, and isovector-tensor densities scaled according 
to NDA: 
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where N is the nucleon field. The lagrangian is also organized according to an expansion in 
derivatives acting on these densities; NDA dictates that each derivative is scaled by A: 

P = ^ . (6) 

Except for the kinetic term, derivatives acting on an individual nucleon field, rather than 
on a density, are eliminated by field redefinitions as in Ref. in favor of terms of the form 
of Eqs. d5)~@- This is because time derivatives on the nearly on-shell valence nucleons are 
of 0(M) and are therefore not suppressed according to (^). As discussed in Ref. 0, this 
procedure is only able to eliminate derivatives in the combination 7 M <9 M (and therefore d^d^). 
We can transform away mixed-derivative terms such as (Nd^N)(Nd^N) in favor of tensor 
terms, which are then neglected (see below) ||. 

Naive dimensional analysis provides an organizational principle that directly translates 
into numerical estimates at the mean-field level. For example, each additional power of 
p s is accompanied by a factor of fA- The ratios of scalar and vector densities to this 
factor at nuclear matter equilibrium density are between 1/4 and 1/7 [0, which serves as 
an expansion parameter. Similarly, one can anticipate good convergence for gradients of 
the densities, since the relevant scale for derivatives in finite nuclei should be roughly the 
nuclear surface thickness a, and so the dimensionless expansion parameter is 1/Acx < 1/5. 
This expansion is only useful, however, if the coefficients are not too large. 

In effective lagrangians of QCD applied to scattering problems, fits to experimental data 
suggest that when NDA is applied the remaining dimensionless coefficients are of order unity. 
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This is known as "naturalness" ; it is an essential feature if NDA is to be useful as an orga- 
nizational scheme. The premise of naturalness in point-coupling models is to be tested here 
for finite-density applications through fits to experimental nuclear data. Naive dimensional 
analysis applied to the point- coupling model of Nikolaus et al. fl8l , |T9H and to the meson 
model of Ref. |9j already suggest that effective nuclear models are natural, which in turn 
implies a convergent mean-field expansion in density based on an organization prescribed 
by NDA ||. This is a nontrivial result, because the naturalness assumption implies that 
all the short- distance physics (with scale 1/A) is incorporated into the coefficients of the 
effective lagrangian, while long-distance finite-density effects should be calculated explicitly 
in a systematic application of the effective lagrangian. At the mean-field level, however, we 
also approximately absorb long-distance many-body effects from ladder and ring diagrams 

The effective lagrangian should in principle contain every possible term (allowed by sym- 
metries) to a given order under this organization. However, certain classes of terms will be 
poorly determined by fits to bulk nuclear observables. Here meson-exchange phenomenology 
and experience with relativistic mean-field meson models are useful guides. Thus, we follow 
the physics motivation of Ref. || to determine which terms will be essential and which can 
be omitted. In particular, each term we include corresponds to one in the meson-nucleon 
lagrangian, as identified through a simple leading-order analysis. For example, at leading 
order, the scalar field, 0, in the meson model of Ref. || is proportional to the scalar density, 
p s ; terms second order in the scalar density here (including those containing derivatives) 
are therefore related to the mass and kinetic terms of the scalar meson field as well as its 
Yukawa coupling to the nucleon in the meson model. The term cubic in the scalar density 
in the point-coupling model has a correspondence with the term cubic in the scalar field, 
and so on. The absence of a tensor boson in the mesonic model corresponds to the absence 
of terms of the form [Na^N] 2 and [daNa^N] 2 in our point-coupling lagrangian. Tensor 
meson masses are large and tensor mean-field densities are small, so even if the coefficients 
are natural, we expect that such terms would have a small effect. An analysis of terms not 
included here is postponed to a future investigation. 

The resulting point-coupling lagrangian is divided into four parts: 



£ — £nn + £tt + £em + ^VMD • 
The pure nucleon contact interactions are contained in £nn, and take the form 
£nn = N[id^ — M)N 



(7) 



+ [d^f v } 2 [( d + a 2 p i 



(8) 



The terms are organized to manifest the expansion and truncation of £ in powers and 
derivatives of the densities. The point-coupling parameters have a "tilde" over them to 
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distinguish them from the corresponding parameters of the hadron lagrangian in Ref. ||. 
We follow Ref. ll| by not including counting factors in applying NDA, although they were 



included in Ref. ||. This prescription will be tested empirically by our fits. 

Compared to the point-coupling lagrangian of Nikolaus et al. Ul8| , our lagrangian includes 



additional terms with coefficients 771, fj 2 , 5i, a.2, f p , and f v . We also exclude contributions 
that would correspond to an isovector, scalar channel in the meson lagrangian. This meson 
was not included in Ref. || based on meson-exchange phenomenology^ and we note that the 
corresponding point-coupling coefficient in [18[] was found to be unnaturally small. Finally, 



we also postpone consideration of terms with an explicit dependence on the four-velocity 
of the nuclear medium (such as u^j!^), which might arise in an effective mean-field energy 
functional. 

The electromagnetic kinetic and interaction terms are contained in £em and £vmd- 
Electromagnetic observables are calculated as an expansion in the electric charge e, as well 
as a derivative expansion. Here we work to first order in e; thus, only couplings to the 
nucleon that are linear in the photon field are considered. The lowest order terms in a 
derivative expansion are contained in £em and take the same form as the photon-nucleon 
coupling terms in the meson model of Ref. Q: 

Cem = -pfA» llM (l + r 3 )N - -JLf^NXoVN 

-^d v F^N{0 s + fashjN - -F^F, V . (9) 
The anomalous magnetic moment, A, is given by 

X = h P (l + r 3 ) + h n (l-r 3 ) , (10) 

where A p and A n are the anomalous moments for the proton and neutron. 

In a mesonic model, an adequate description of the low-momentum electromagnetic form 
factors of a single nucleon is achieved through a combination of vector-meson dominance 
(VMD) and direct couplings of the photon to nucleons PP3|. Since heavy-meson 

degrees of freedom are eliminated here in favor of contact interactions, these effects must 
be expressed in the point-coupling model through direct higher-derivative couplings of the 
photon to the nucleon, contained in £vmd- I n practice we only include the next correction 
to the charge form factors: 

£vmd = -^d 2 d u F^N({5 s + S v r 3 } lfl )N . (11) 

The coefficients of Cem and £vmd are fixed by experimental electromagnetic data for protons 
and neutrons in the vacuum and are not varied independently when fitting to nuclear data. 
This is detailed below. The end result is that the composite electromagnetic structure of the 



1 The NN interaction in that channel is weak plf . There is no meson with these quantum numbers 
with a mass below 1 GeV, and two (identical) pions in a J = state cannot have T = 1, so there 
is no analog to the a meson. 
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nucleons is incorporated as a derivative expansion, which is appropriate for low-momentum 
physics. In the point-coupling model the expansion is constructed explicitly (at tree level) 
to a given order in momentum. In a meson model, vector meson dominance incorporates 
contributions to the form factor to all orders in the momentum. 

The pion kinetic and interaction terms are contained in L^. The procedure for construct- 
ing this part of the lagrangian is similar to that presented for the meson model §J. Since 
the pion field vanishes in the mean-field approximation, the details of its construction are 
not presented here, nor is an explicit form for these terms given. The pion will first enter 
when considering two-loop contributions to the energy. 

Beyond a leading-order transformation from a meson lagrangian to a point-coupling 
lagrangian, the full expressions for the mean meson fields involve infinite series in powers of 
the various bilinears of the nucleon field. A precise transformation from the meson effective 
lagrangian at the mean-field level would therefore lead to the presence of higher-order terms 
in a point-coupling model that are treated here as negligible. In practice, because of delicate 
cancellations, a simple truncation of these terms leads to a poor fit with experimental data, 
but if the parameters are allowed to readjust slightly the new fit to the data is quite good. 

This is not unlike the truncation of higher-order terms within the mean-field meson 
model: the inclusion of fifth-order terms (see Eq. (54) of Ref. ||) into the lagrangian im- 
proves a fit to the data only marginally, but if the fifth-order terms are then simply truncated 
from the lagrangian, their effects must be absorbed into an adjustment of the lower-order 
coefficients if a good fit is to remain intact. Thus, the correspondence between the point- 
coupling model used here and the mean-field meson model of Ref. || is not exact due to 
truncations. Relations of the leading-order point- coupling coefficients here to the corre- 
sponding mean-field meson model coefficients that would result from a full transformation 
are listed below in Tables |TJ, |TTT|, [TV| and 0. 

In the meson model of Furnstahl et al. 0, the vector and p-meson masses are fixed at 
their vacuum values of 782 MeV and 770 MeV, respectively. The correspondences of these 
masses to the point-coupling parameters are given by 

ml = ( 2 /( d , m 2 p = l 2 /l d . (12) 

The vector meson and p-meson masses can therefore be "held fixed" within the point- 
coupling model by varying only ( 2 and £2 and determining ( d and ^ through these relations. 
These relations were not imposed in the point-coupling model studied by Hoch et al. [JTS 



Here we perform fits to experimental data both with these combinations fixed and allowing 
the four parameters to vary independently. 



III. SINGLE-NUCLEON PROPERTIES 

The values of the photon coupling parameters A n , X p , (3 S , /3 V and 5s+^v are all determined 
from single-nucleon electromagnetic form factors in the vacuum, calculated at tree level. The 
anomalous magnetic moments of the nucleon are fixed at X p = 1.793 and A n = —1.913 0. 
The other three parameters are related to the isoscalar and isovector charge form factors, 
which are given here for spacelike momenta Q 2 = — q 2 by 
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F '^) = \-m-^-w^ + -- < 14) 

The coefficients of the second-order terms are proportional to the corresponding mean-square 
charge radii, and their values are fixed from the experimental values 

f3 s = M 2 (r 2 ) sl /6 « M 2 (0.79fm) 2 /6 = 2.36 , (15) 
P v = M 2 (r 2 ) vl /Q « M 2 (0.79fm) 2 /6 = 2.36 . (16) 

These coefficients combine the direct coupling and VMD contributions that together de- 
termine the charge radii in the model of Ref. ||. Unlike that model, there is no density 
dependence here from the effective masses of the vector mesons in medium. Since the dif- 
ference between the point proton density and the full charge density is not very important 
in determining the self-consistent wavefunctions and energy levels, the point-coupling elec- 
tromagnetic contributions produce similar results to the conventional convolution procedure 
|nj, in which empirical single- nucleon form factors are folded with point nucleon densities 
after self-consistency is reached. 

Nuclear structure observables do not depend strongly on the values of the fourth-order 
coefficients S s and 5 V . Based on the criteria used here, the quality of fits are not affected by 
the inclusion of these parameters and their associated terms. They do play a key role in the 
(momentum-space) charge form-factor of each nucleus for higher values of Q 2 , however. The 
d.m.s. charge radius, which depends on the position of the first zero of the charge form-factor 

R dms = 4.493/QW , (17) 

is used as one criterion for optimization, and the parameter sets with good values of \ 2 
do reproduce the experimental result for the value of Qq . However, if the fourth-order 
corrections are omitted, for oxygen there is a significant deviation from experiment beyond 
this value of momentum. In particular, the second maximum falls short of the predicted value 
(see Fig. p. This is in contrast to the meson model, where vector-meson dominance yields 
the necessary momentum dependence to produce accurate results for the second maximum 
without including it as one of the criteria for optimization. 

Fourth-order corrections in momentum are therefore included in the model through the 
parameters S s and 8 V . A determination of the values of these parameters through vacuum 
properties is sufficient to reproduce the second maximum (see Fig. [I]). In fact, due to the 
small magnitude of the neutron charge density in comparison to the proton, only the isospin 
+ 1/2 components of the additional terms are needed and 5 S and 5 V are therefore set equal. 

We use a dipole fit to the proton form factor to determine the value of 5 S = 8 V . (Any 
alternative parameterization of the form factor could be used instead; our results are not 
sensitive to the details because it is a low- momentum expansion.) The proton form factor, 
G P E (Q 2 ) is known to be fitted well by a dipole form [25|: 



G P E (Q 2 )=[l + *^) =1-2^ + 3^ + .... (18) 
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FIG. 1. Charge form factor of 16 O. The solid line is taken from Ref. |24| . Form factors are shown 
for set FZ3 of the point-coupling model with and without fourth-order momentum corrections. 
Also, set Gl from the meson model Pi is shown. 



This form factor is related to the isoscalar and isovector form factors by the relation 



G p E {Q 2 ) = Ff (Q 2 



AM 2 



F$(Q 2 ) 



where 



f?(q 2 ) = f:(q 2 ) + f»(q' 



(19) 



(20) 



Using Eqs. (O) and ( [Hp and the fact that 5 S = 5 V , the proton form factor in Eq. flTP| ) can 
be expanded as 



G p E (Q 2 



1 - 



2f3 s + 2f3 v + Xp 2 _ ^s_ n 4 , 
W M 4 ^ 



Comparison with Eq.(ffH) yields 



a 



4M 2 



2& + 2& + A p 



(21) 



1.39 ; 



28, = -3a 2 = -5.80 



(22) 
(23) 



In principle, higher-order momentum corrections to the anomalous form factors F£ ,v , 
which are not considered here, would lead to an additional Q A dependence to the proton 
form-factor. Our sole interest is in determining the overall coefficient appearing at fourth- 
order in Eq. (|21|), however, and such a contribution could be re-absorbed into a redefinition 
of <L 
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IV. THE MEAN-FIELD SOLUTIONS 



We work at one-baryon-loop order in this paper, which is equivalent to the Dirac-Hartree 
approximation ||10|| . This should be sufficient to test the consistency of naturalness through 



NDA and the truncation of the lagrangian. As discussed in Ref. ||26|| , a Hartree calculation 
can be viewed as equivalent to a density-functional approach, in which higher-order many- 
body corrections are treated approximately. Such corrections have not yet been explored in 
a relativistic point-coupling model. In future work, the stability of the Hartree results will 
be tested in a two-loop calculation. 

From the lagrangian, one can derive the Dirac equation and the energy functional for 
nuclei conventionally, keeping in mind that there are time derivatives acting on ip^ as well 
as ip. (See Ref. [JH] for details.) Because of its large mass, loop integrals with the nucleon 
include dynamics from distance scales that are much shorter than the scale set by the valence- 
nucleon momenta, which are limited by the Fermi momentum k-p- These short-range effects 
are included implicitly in the coefficients of the lagrangian and energy functional. Formally, 
one can include counterterms to remove these loop effects to all orders, which is always 
possible, since all such terms are already contained in the effective lagrangian. 

The single-particle Dirac hamiltonian for spherically symmetric nuclei takes a form sim- 
ilar to that of Ref. 0: 

h{x) = -ia-V + W(r) + ^r 3 R(r) + (3{M - $(r)) + i(3a ■ T(r) + l -(3a ■ T 3 (r)r 3 
+\0- + rs)A(r) + ^(A + (3 v r 3 )V 2 A(r) - -^(5 S + ^r 3 )(V 2 ) 2 A(r) 

~iu [Xp{1 + Ts) + An(1 " T ^ a ' VA ^ > ( 24 ) 

where the potentials $, W, R, A, T, and T 3 are defined by 

$ = -A{2R 2 p s + 2n d V 2 p s + 3£ 3 p| + 4^ + Ujv) 2 + ^hQvfh + V P Gl? 

+ a^Vpsf + 2a 1 p s V 2 p s - 5 2 (Vj°) 2 } , (25) 
W = A{2C 2 (j°) + 20 V 2 (J°) + 4C 4 (j°) 3 + 2^J°p s + 2fj 2 ]° v p 2 s 

+ 2a 2 (Vp s -V]° v + p s V 2 J°) + f v V-s} , (26) 

R ee A{2&3 + 2i d V 2 ]° T + 2rj p p s ]° T + f p V-s T } , (27) 

A = eA ^ (28) 

T ee -A/ B V3S , (29) 

T 3 ee -Af p Vf T , (30) 

Si EE S 0i , (31) 
s -ri = S-rOi ■ (32) 

The Dirac equation with eigenvalues E a and eigenfunctions A Q (x) is 

hN a (x) = E a N a (x) , J d 3 x At(x)A Q (x) = 1 . (33) 
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Following the conventions of Ref. |1C| (also used in Ref. 0), the eigenfunctions for spherically 
symmetric nuclei are written in terms of spin spherical harmonics, $ Km : 

"H^!:!t>* (34) 

where t = 1/2 for protons and t = —1/2 for neutrons. The equation for N a then reduces to 
a set of coupled equations for G a and F a : 

4- + - J G a (r) - [E a - U^r) + U 2 {r)]F a {r) - U 3 (r)G a (r) = , (35) 



dr r 
' d k 
dr r 



+ [E a - U^r) - U 2 (r)]G a (r) + U 3 (r)F a (r) = , (36) 



where we have defined the single-particle potentials by 

U x (r) = W(r) + t a R(r) + (t a + l/2)A{r) + ^ (ft + 2t a p v )V 2 A(r) 

' (5 s + 2tJ v )(V 2 ) 2 A(r) , (37) 



M 4 

U 2 (r) = M- $(r) , (38) 

1 dA 
U 3 (r) = -T(r) • f - f B T,(r) • f + —[(A, + X n )/2 + t a (X p - X n )} — . (39) 

The mean-field densities can be expressed in terms of the radial wave functions G and 
F in the same manner as Refs. |J and [TO] : 



(Gl{r) ~ F 2 a (r)) , (40) 



m = ^-^(G 2 a (r) + F 2 a (r)), (41) 

S(0 = ^S^ 1 (O(^.(r) + ^(r)) ) (42) 
1 occ 2i + 1 

^•^^E^T^GW), (43) 

«r(r) - = «T E 2G a (r)F a (r) . (44) 

The summation superscript "occ" means that the sum runs only over occupied (valence) 
states in the Fermi sea. The mean-field equation for the photon field is given by 

- V 2 A = e 2 Pchg (r) , (45) 

where the charge density is 

OCC 1 1 ~ 

~ ^i(V 2 ) 2 iV Q (^ + 5 v r 3 ) l0 N a + ^ V ■ N a X(3aN a } . (46) 
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The Hartree equations are solved by an iterative procedure similar to that used for 



meson models ||10|| . The procedure is simplified over that of the meson model in that only 
one non-linear differential equation for a meson field (the Coulomb field A) is solved; the 
input potentials for each iteration of the Dirac equation are then just evaluated in terms of 
the densities from the previous iteration. See Ref. [|18j for more details. 

An expression for the energy is given by the first diagonal element Too of the mean-field 
stress-energy tensor: 

E = Jd 3 xT 00 

/OCC n S 

d 3 x N a {-i(3a.-V + M)N a + / 2 A 2 J d 3 x ^ 2 p 2 s - K d [(d°p s ) 2 + (Vp s ) 2 ] 

+ K 3 p 3 s + £ 4 Ps + ViCfvfps + rhGv) 2 ?s 



+ Uf T f - a(£>°j£) 2 + (vj?) 2 ] + v P &)%_ 

- a.pslidops) 2 + (Vp s ) 2 ] - a 2 p s {0o] O v ) 2 + (Vj°) 2 ] 
+ f v V~f v -s + f p V] T -8^ 

+ electromagnetic terms . (47) 

Ground-state densities are time- independent, so all time derivatives in E vanish. The Dirac 
equation can be used to replace the first two terms with a new expression involving a sum 
over the energies of the occupied states; after integrating by parts and using the equation 
for the electromagnetic potential, the expression becomes 

E = Y,E a - d 3 x f 2 A 2 l K 2 p 2 s + K dPs V p s + 2k 3P 3 s + 3n 4 p s 

a ^ 

+ Ufv) 2 + OjvV 2 J° + 3C 4 (J°) 4 
+ 6(l°) 2 + ^v 2 l° + 2ry p (l°) 2 p s 

- 2a 1 p s {Vp s ) 2 - 25 2 p s (Vj°) 2 

+ f p ]° T V • s T + f v j° v V ■ s} - l - Jd 3 x Ap chg . (48) 

A center-of-mass correction for the energy is incorporated as in Ref. but the charge radius 
is not corrected (including the latter does not change our conclusions). 



V. OPTIMIZATION 

All parameters of the model other than those fixed by single-nucleon properties in the 
vacuum (see Sect. Ill) are treated as free variables and an optimum fit to experimental 
nuclear data is sought. The same observables and weights are used as in Ref. 0. This 
optimization faces serious challenges. The existence of delicate cancellations at lowest order 
in the density expansion necessitates calculations to higher order in the density to obtain an 
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adequate description of nuclei. This increases the size of the parameter set to optimize and 
since the observables themselves are highly correlated, there is a problem with underdeter- 
mination. The delicate cancellations also result in an extreme sensitivity to the lower-order 
parameters, precluding large steps through the parameter space if the Hartree iterations are 
to remain stable; the optimization program must therefore be efficient at navigating narrow 
valleys. 

Nonetheless, the optimization procedure has been significantly improved over that used 
in obtaining the results in Ref . || . This has been achieved in part by exploiting the relatively 
small difference between scalar and vector densities at ordinary nuclear densities. We can 
rewrite the lagrangian as an expansion in the bilinears p + and p_ defined byQ 

p+ = L = ^(ps + Jv); (49) 

„ _^ N a (l- lQ )N a 1 ^ 

P-=2^ =2^ ps ~ 3 ^' ( 5 °) 

where the sum is over occupied states N a . This separation is reminiscent of the heavy baryon 
formalism (HBF) |f27| , |28| in the vacuum, with the spin matrix (1 — 7o)/2 acting to project 
out the negative energy states in the vacuum (at least to leading order in an inverse nucleon 
mass expansion). 

Here we simply observe that a useful hierarchy for optimization is given in terms of 
"optimal parameters," which are listed in Table |. These parameters are linear combinations 
of the couplings in the lagrangian. The hierarchy is based on the observation that the 
difference p_ between scalar and vector densities is small and empirically scales like p^/ 3 near 
equilibrium density^ We can classify terms according to a power of p + , which is roughly 
1/4 to 1/7. The label in the first column of Table p] indicates how different parameter sets 
are organized (see below); it does not correspond precisely to a systematic expansion in p + . 
Performing the optimization in terms of the fli parameters in Table [I] is more efficient; we 
stress, however, that the energy functional has not been changed. 

Improvements were also made to the optimization procedure by keeping careful track of 
the necessary precision of the parameters. The calculation of x 2 (treated as a function of the 
optimal parameters) is much more sensitive to small changes in the lower-order parameters 
than changes of the same size in the higher-order parameters. The difference is severe 
enough to create problems for the minimization software. To remedy this, each parameter 
was expressed in the form 

tt t = of* ed + 10-^5Qi , (51) 

where fif xed and A, > are determined by the initial values of the parameters and are 
held fixed during optimization. The parameters are then varied by changing the value of 



2 An analogous procedure can be performed in the meson model M by rewriting the lagrangian in 
terms of the sum and difference of the scalar and (zero component of) the vector fields, (m 2 s Q/gl ± 
mlW/g 2 v )/2. 

At low density p- scales like p! , but this density dependence is independent of the parameters. 
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TABLE I. Hierarchy of "optimal" terms (see text). Terms are presented in order of importance 
from top to bottom. Derivative terms are treated separately. 



Without derivatives 


Terms in lagrangian 


set 


optimal parameter 


new notation 


covariant 


Po 
Po 


1 =K 2 + (2 
0,2 = K 2 - (2 






Pi 

P2 


3 = K 3 + 771 

4 = K 3 - 77i/3 


n 3 (p% + pi) + 3^ 4 (/^p- + p+pI) 


«3£| + ?yi(i°) 2 p s 


P3 
Pa 
P5 


o 5 = «4 + C4 + m 

66 = K A — (a 

67 = ^4 + (a - V2/3 


n 5 (p\ + pi) + 40 6 (p3 p_ + p +/ ?3 ) 

+ 6Q 7 (p 2 + p 2 _) 


kap 4 s + CA(j°v) 4 + mp 2 s Cfv) 2 


With derivatives 




^0 
Po 


Ai = n d + Q 

A 2 = Kd~ Cd 


Aip M p + ) 2 + (^_) 2 ] 


Hdidf.p.s) 2 + (s(d fl jvu) 2 


P2 
P2 


A\ = Oi\ + OL 2 

A 2 = a\ — a 2 


yl 1 P / ,p + ) 2 + (^p_) 2 ](p + + p_) 
+ 2l 2 ^p+^p_(p + + p_) 


on{dp,p s ) 2 Ps + Oi 2 (d^i vv ) 2 p s 



<5f2j. Lower-order parameters have a larger value of \ so that large changes in 5Qi actually 
correspond to small changes in the true parameter The values of fij and Aj were chosen 
to yield roughly similar curvatures of the x 2 ~function for variations of each 5Qi and therefore 
yield a similar sensitivity of x 2 to each parameter. Given the delicate dependence on the 
lower-order parameters, an advantage was also gained over the optimization procedure used 
in the meson model by creating a routine outside the minimization package to calculate 
the gradient of the x 2 ~function (with respect to the optimal parameters) needed by the 
optimizer. 

VI. RESULTS AND DISCUSSION 

The model is fit to experimental data for four different categories of parameter sets, 
each at different levels of truncation. The parameter sets are categorized according to three 
alphanumeric labels. Optimizations in which the fixed-mass relations from the meson model 
are imposed on the parameters (d and ^ are labeled with an initial letter of "F" , and those 
in which Q and are freely varied are labeled with an initial letter of "V". Optimizations 
in which the parameters ol\ and a 2 are fixed at zero are labeled with a second letter of 
"Z" , and those in which they are freely varied are labeled with a second letter of "A" . The 
nomenclature for each optimized parameter set also includes a number reflecting the level 
of truncation (see Table |) : 
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TABLE II. 


Parameter sets from fits to finite nuclei. Fixed 


771,, 771 „ 

' ' V V 3 " "p 


and ol\ = ot 2 


= 0. 


Point-Coupling 


Meson 


FZO 


FZ1 


FZ2 


FZ3 


FZ4 




-g 2 f 2 /(2m 2 ) 


-1.930 


-1.229 


-1.566 


-0.963 


-1.042 






1.490 


0.760 


1.079 


0.460 


0.538 


6 




0.679 


0.524 


0.664 


1.014 


0.992 


m 






0.521 


-1.021 


1.615 


0.658 


m 










-1.178 


-1.207 








0.174 


1.784 


-0.368 


0.608 












-0.196 


-1.530 


Ci 










-0.196 


1.128 


% 








-1.142 


-3.400 


-3.245 


fv 






0.411 


0.434 


0.667 


0.685 


fp 






2.503 


1.904 


1.659 


1.575 




-g 2 f 2 A 2 /(2mi) 


-3.140 


-2.006 


-2.406 


-1.779 


-1.859 


Q 


g 2 J 2 A 2 /(2mi) 


1.445 


0.737 


1.046 


0.446 


0.522 




g 2 p f 2 A 2 /(2mj) 


0.679 


0.524 


0.664 


1.014 


0.992 


K 2 /Kd 


m 2 /A 2 


(0.784) 2 


(0.783) 2 


(0.807) 2 


(0.736) 2 


(0.749) 2 


-K 2 A 2 /(f 2 K d ) 


ff s 7(4vr) 


1.009 


0.640 


0.866 


0.443 


0.497 


C 2 A 2 /(^Cd) 


9 2 J(^) 2 


1.306 


0.666 


0.946 


0.404 


0.471 




9 2 p /(^) 2 


0.577 


0.445 


0.564 


0.862 


0.843 


x 2 




1600 


65 


41 


37 


37 



parameter set of terms kept 
P = {f2i,J\,Ai,A 2 ,e2,fd} 
Pi = PoU{fi 3 ,/,,/v} 

p 2 = p 1 u{n 4 ,fj p } 

P 3 = P 2 U {Q 5 } 

p 4 = p 3 u {n 6 } 

We stress that the number denoting the level of truncation does not correspond with an 
expansion to that order in the scalar and vector densities, but instead to an organization 
according to the optimal parameters. However, we can identify set Po with a truncation at 
second order in the densities, set P 2 with third order in the densities, and set P4 with fourth 
order in the densities. Sets FAi and VAi also include the parameters Si and a 2 \ only values 
of i > 2 are considered in these categories. 

The optimizations were carried out with the Minuit software package |?9| . The resulting 
data sets for fits at various orders of truncation are presented in Tables [TT], [TT1], [IV] and 
[V[ Results from parameter sets FZO and VZ0 show that truncation at second order in 
density leads to a poor fit with experimental data and yields a similar x 2 to the analogous 
parameter set Wl in the meson model from Ref. 0. A significant difference from the 
meson model arises at the next level of truncation, however: the inclusion of only a single 
optimal parameter beyond second order in density yields a significantly better fit than the 
full third-order results of the meson model 0. Incorporating all third-order terms, sets FZ2 
and VZ2 are comparable to or better than the best fourth-order results of the meson model 
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TABLE III. Parameter sets from fits to finite nuclei. Freely varied m v , m p , ol\ = «2 = 0. 



Point-Coupling 


Meson 


vzu 


V LI 


V LI 


V L6 


^ 7*7 A 




- 5 2 / 2 /(2m 2 ) 


—1.859 


— 1.238 


—1.972 


— 1.509 


— 1.482 


£2 


„2 f2 l/r lrm 2\ 


1.427 


U.7o5 


1 A OO 

1.482 




0.961 




„2 t2 //'o™2\ 


0.596 


n trot' 

0.535 




0.345 


0.521 






U.OZD 


— l.yUo 


l.Uo < 


u.uoo 


m 










— i.yzy 


1 onn 

— i.yuu 








0.175 


2.509 


0.202 


1.278 












-0.321 


-2.105 


a 










-0.321 


1.471 










1.677 


1.133 


-0.220 


fv 






0.412 


0.027 


0.041 


0.189 


fp 






2.836 


2.889 


3.780 


3.267 


Kd 




-0.732 


-2.531 


-1.740 


-0.639 


-0.712 


Cd 


glfl^/(2mi) 


-0.739 


1.238 


0.389 


-0.673 


-0.603 


id 


^/2 A 2 /( 2 m 4) 


-3.813 


1.074 


-1.578 


-3.060 


-2.410 


K2/K d 


m|/A 2 


(1.594) 2 


(0.699) 2 


(1.065) 2 


(1.537) 2 


(1.443) 2 


C2/Cd 


m 2 /A 2 


-1.931 


0.620 


3.811 


-1.468 


-1.593 




m 2 /A 2 
5 2 /(16vr 2 ) 


-0.156 


0.498 


-0.164 


-0.112 


-0.216 


-2~4^/utK d ) 


4.014 


0.515 


1.900 


3.031 


2.625 




5,V(16vr 2 ) 


-2.342 


0.404 


4.801 


-1.233 


-1.300 


^W/uUd) 


^/(16vr 2 ) 


-0.079 


0.226 


-0.036 


-0.033 


-0.096 


x 2 




1230 


62 


35 


27 


26 
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TABLE IV. Parameter sets from fits to finite nuclei. Fixed m v , m p . Nonzero a\ and a 2 . 



T"» * j_ /"I 1 ' 

Point-Coupling 


Meson 


1A2 


1A3 


T71 A A 

1A4 




-<7 2 / 2 /(2m 2 ) 


—2.065 


— 1.306 


— 1.298 


> 

C2 


2 -f2 llr, 2\ 


1.563 


r\ '~7or\ 

0.789 


0.782 




g 2 p f 2 /(2m 2 ) 


0.466 


1.050 


1.045 


m 




no 

—2.512 


1.007 


0.665 


m 






— 1.513 


— 1.547 


^3 




3.200 


0.276 


0.625 








-0.252 


-1.125 


a 






-0.252 


0.609 






0.107 


-3.926 


-3.864 


/„ 




0.000 


0.315 


0.342 


/p 




2.449 


1.857 


1.856 


Kd 


-g 2 f 2 A 2 /(2mf) 


-3.117 


-2.190 


-2.136 


Cd 


g 2 f 2 A 2 /(2mt) 


1.515 


0.765 


0.758 


id 


g 2 f 2 A 2 /(2m*) 


0.466 


1.050 


1.045 


OL\ 




16.42 


11.138 


9.125 


OL2 




-13.12 


-9.818 


-8.387 




m 2 s /A 2 
5 s V(16vr 2 ) 


(0.814) 2 


(0.772) 2 


(0.779) 2 


-K 2 2 /(87T 2 K d ) 


1.162 


0.662 


0.670 


C2/(8n 2 Q) 


5, 2 /(16vr 2 ) 


1.372 


0.692 


0.686 


tt/^Hd) 


5 2 /(16vr 2 ) 


0.396 


0.893 


0.888 


x 2 




36 


33 


32 
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TABLE V. Parameter sets from fits to finite nuclei. Freely varied m v , m p , a.\ and «2- 



Point-Coupling 


Meson 


VA2 


VA3 


VA4 


K 2 


-3 2 / 2 /(2m 2 ) 


-2.064 


-1.471 


-1.192 




o 2 / 2 /(2m 2 ) 


1.563 


0.948 


0.667 




g 2 J 2 /(2m 2 ) 


0.263 


-0.048 


0.364 


/ -L 




-2.519 


1.172 


1.987 








-1.885 


-2.445 


^3 




3.209 


0.114 


-0.501 


K4 






-0.314 


-1.607 


C4 






-0.314 


0.792 


rj p 




1.670 


4.597 


1.148 


fv 




0.023 


0.142 


0.251 


fp 

■J fj 




2.888 


3.552 


3.328 


Kd 


-g 2 f 2 A 2 /(2mi) 


-3.066 


5.019 


4.675 


Id 


5 2 / 2 A 2 /(2m4) 


1.474 


-6.310 


-5.812 




g 2 f 2 A 2 /(2m 4 p ) 


-1.456 


-4.954 


-3.632 


Oil 




15.63 


-37.83 


-37.43 


5(2 




-12.62 


38.14 


35.84 




m 2 s /A 2 


0.673 


-0.293 


-0.255 




ml/A 2 


1.060 


-0.150 


-0.115 




m 2 /A 2 


-0.180 


-0.009 


-0.100 


-Ki/(87T 2 K d ) 


9 2 /(16^ 2 ) 


1.182 


-0.367 


-0.258 


C 2 /(87T 2 0) 


5, 2 /(16^ 2 ) 


1.408 


-0.121 


-0.065 


£/(87r 2 &) 


5 2 /(16^ 2 ) 


-0.040 


-0.0004 


-0.031 


x 2 




33 


25 


23 



TABLE VI. Optimal coefficients for selected parameter sets. 



parameter 


order 


FZ2 


VZ2 


FA2 


VA2 


FZ4 


VZ4 


FA4 


VA4 


Hi 


f + 


-0.49 


-0.49 


-0.50 


-0.50 


-0.50 


-0.52 


-0.51 


-0.52 


^3 


P% 


0.76 


0.60 


0.69 


0.69 


1.27 


1.33 


1.29 


1.49 


^2 


P+P- 


-2.65 


-3.46 


-3.63 


-3.63 


-1.58 


-2.45 


-2.08 


-1.86 


5 


P\ 


0.00 


0.00 


0.00 


0.00 


-1.61 


-2.53 


-2.06 


-3.26 




f+P- 


2.12 


3.14 


4.04 


4.05 


0.39 


1.26 


0.40 


-1.16 


^6 


P+P- 


0.00 


0.00 


0.00 


0.00 


-2.66 


-3.58 


-1.74 


-2.40 


VL 7 


pIp 1 - 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


Ai 


Pl 


-1.36 


-1.35 


-1.60 


-1.59 


-1.34 


-1.31 


-1.38 


-1.14 


A 2 


p+p- 


-3.46 


-2.13 


-4.63 


-4.54 


-2.38 


-0.11 


-2.90 


10.49 


A 1 


p\ 


0.00 


0.00 


3.30 


3.01 


0.00 


0.00 


0.74 


-1.59 


A 2 


~2 ~ 

p+p- 


0.00 


0.00 


29.54 


28.25 


0.00 


0.00 


17.51 


-73.30 
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FIG. 2. Best-fit x 2 values at different levels of truncation (see text). Results for the lowest 

level of truncation are not shown. 



given in Ref. ||. Because of the improvements in the optimization procedure, however, a 
quantitative comparison of x 2 values from this work and from the meson-model calculations 
of Ref. [EJ might be misleading. In a future investigation the meson model will be reanalyzed 
to determine if the differences noted here are due to the improvements in optimization or 
arise from a greater flexibility of the point-coupling model in fitting nuclear properties. 

Beyond the lowest level of truncation, fits to the charge radii, d.m.s. radii, binding 
energies, and spin-orbit splittings for each data set are all quite good; they are essentially 
all at the relative accuracy prescribed by the corresponding weights (see Ref. || for details). 
The fits here with i > 2 are as good as or better than the best fits in Ref. || . Also, the low- 
momentum behavior of the nuclear charge form factors is well reproduced and is comparable 
to those of the meson model |J. Since our focus here is on testing naturalness and the NDA 
expansion in the point-coupling model, we leave a detailed comparison between meson and 
point-coupling models to a future work. 

A plot of the value of x 2 f° r parameter sets in each category at different levels of trunca- 
tion is shown in Fig. |2|. The difference in x 2 between the third and fourth level of truncation 
is small and is a signal that including higher-order terms will not improve the fit to data 
significantly. This has been verified by optimizing fits for parameter sets containing the 
leading fifth-order optimal parameter. Also, at each level of truncation, parameter set VZ 
yields a better fit than FA, and FZ and VA yield the worst and best fits, respectively. That 
is, removing constraints on the parameters systematically improves the fits. Note however 
that all of the fits for i > are good! We should be cautious, therefore, about drawing 
strong conclusions based on such small differences in the \ 2 ■ 
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A. Naturalness 



A review of the parameter sets FZi and VZz shows that most of the coefficients are 
natural. For sets in which Si and fi 2 are non-zero (FAi and VAi), however, these parameters 
tend to be large and unnatural. This tendency can be understood by considering the optimal 
parameters defined earlier. They are given in Table |V| for sets with i = 2 and 2 = 4. The 
second column indicates the leading power of p + and p_ for the corresponding term. 

The coefficients are ordered according to the empirical rule that p_ should be counted 
roughly as fff 3 (near equilibrium density). We expect that lower order parameters will 
be better determined by the fit to the data, and this is indeed the case. (In general, the 
optimal parameters are much better determined than the original covariant parameters.) 
Furthermore, when a set is truncated, the highest-order nonzero coefficients are forced to 
compensate for the omitted terms. In sets with 2 = 3, the values of Clx, f2 2 , ^3 and are 
consistent with the same range in values as the sets with i = 4. Truncation of the term 
corresponding to Qq can be compensated by Q§ (which is underdetermined at i = 4 anyway). 

The truncation at i = 2 is different, however. These sets give quantitatively good fits to 
the data, but although the fits yield consistent values of fi 2 and Cl 3 for these sets, their values 
are significantly different from those in sets with i = 3 and i = 4. In fact, the large size 
of f2 2 compared to Q% for i = 2 indicates a breakdown of the systematics: the contribution 
to the binding energy from the term corresponding to Jl 2 is larger than for f^, while the 
opposite behavior is expected. Thus, although the i = 2 level of truncation yields good fits, 
it is clear that a higher level of truncation is necessary for consistency and to produce the 
systematics implied by naturalness. 

From the table, it is clear that the optimal parameter A 2 is unnatural, while the optimal 
parameter A\ = a,\ + <5 2 is of order unity in magnitude. A similar situation exists for A 2 = 
2(k^ — Cd) and Ai = Kd + Cd'- whereas the parameter Ai seems to maintain a consistent value 
on the order of —1.3 throughout the sets, the value of A 2 varies wildly. The combination 
Ai + Ai(p s ), where (p s ) is an average value (roughly 1/10) is even more tightly determined. 
Since A 2 is poorly determined by the observables, it is not surprising that a fit could lead 
to individual values for a\ and a 2 that are unnatural. 

The point-coupling model is evaluated here for naturalness in the absence of any combi- 
natorial counting factors (as in Ref. [|18j). This is in contrast to the meson model analysis, 
where a term with n powers of the scalar or vector field incorporated factors of 1/n! in the 
treatment of naturalness ||. Yet each analysis yielded natural coefficients in essentially all 
cases. We have not yet understood this difference; a reanalysis of the meson models using 
the improvements to the minimization procedure discussed here is needed before drawing 
firm conclusions. 



B. Nuclear Matter 

Saturation properties of nuclear matter for each of the parameter sets are given in Ta- 



ble [VI]] . The saturation point and binding energy per nucleon have values characteristic of 
relativistic mean-field meson models that successfully reproduce bulk properties of nuclei 
The consistency of the numbers for these sets and the ones in Ref. [HI is striking. Dif- 



ferent weights might shift the numbers slightly, but within small errors the binding energy 
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TABLE VII. Nuclear matter equilibrium properties 



Set 


E/B - M 




K 


04 


M*/M 


Vector Potential/M 


FZO 


-16.1 


1.27 


570 


39.4 


0.533 


0.387 


FZ1 


-16.0 


1.30 


390 


32.6 


0.697 


0.234 


FZ2 


-15.9 


1.30 


320 


31.8 


0.676 


0.253 


FZ3 


-15.9 


1.31 


300 


29.9 


0.742 


0.191 


FZ4 


-15.9 


1.31 


300 


30.0 


0.743 


0.190 


FA2 


-16.0 


1.30 


290 


33.5 


0.601 


0.322 


FA3 


-16.0 


1.31 


280 


29.7 


0.677 


0.252 


FA4 


-16.0 


1.31 


290 


29.7 


0.681 


0.248 


vzo 


-16.3 


1.29 


560 


37.9 


0.535 


0.383 


VZ1 


-16.0 


1.30 


380 


32.8 


0.698 


0.234 


V ZJZj 


— 16 


1 30 


330 


34 8 


597 

U . Oct 1 


397 


VZ3 


-16.0 


1.30 


300 


34.6 


0.623 


0.303 


VZ4 


-16.0 


1.31 


310 


32.8 


0.645 


0.282 


VA2 


-16.0 


1.30 


290 


35.0 


0.602 


0.322 


VA3 


-16.1 


1.30 


300 


39.3 


0.629 


0.297 


VA4 


-16.1 


1.30 


320 


34.5 


0.659 


0.269 



per nucleoli is 16MeV and the Fermi momentum is 1.3 fm 1 . The compressibilities are of 
order 300 MeV for the more complete models, which is somewhat high compared to the best 
fits from Ref. ||. 

The scalar and vector potentials are in many cases much smaller than typical values 
found in meson models | 26fl . Previous experience with meson models implied that large 



potentials were needed to reproduce experimental spin-orbit splittings. The point-coupling 
model of Ref. [18], with M* /M of 0.58, was consistent with this result. 

The new feature here that allows for this deviation from conventional wisdom is the 
isoscalar vector-tensor coupling proportional to /,,, the analog of which has not usually been 
included in relativistic mean-field modelsf] (often based on arguments from vector dominance 
and the smallness of the isoscalar anomalous moment of the nucleon). It provides, however, 
an independent contribution to the spin-orbit force. In Ref. |§, where the analog of this 
term (proportional to f v ) was included, the optimization did not lead to particularly large 
values of the tensor coupling /„, although for set G2 it was large enough to drive the value of 
M* /M to 0.66. This value is outside the range 0.58 < M*/M < 0.64 for conventional models 
with good fits to nuclei p6(] . In the point-coupling model there is apparently a significant 
advantage to having larger /„ (at least with our choice of observables and weights) and we 
find M*/M as large as 0.74! 

Since we are primarily concerned here with testing naturalness and the NDA expansion 
and not with generating parameter sets for practical use, we have simply checked that 



4 We note that the parameter f v in the meson model plays a role in determining electromagnetic 
properties through VMD effects, whereas f v in the point-coupling model does not have any direct 
association with such effects. 
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FIG. 3. Nuclear saturation curves for point-coupling models VA4 and FZ4 from this work, the 
point-coupling model from Nikolaus et al. [18|, and the meson model Gl from Ref. ||. 



setting f v = and re-optimizing the parameters yields a smaller value for M* (consistent 
with conventional wisdom of nuclear models) and does not change any of our conclusions. 
This is indeed the case. A more complete discussion of the spin-orbit force and tensor 
couplings in chiral effective field theories will be given elsewhere |30| . 



In Fig. RL the saturation curves for sets FZ4 and VA4 are plotted against the results of set 



Gl from the meson model of Ref. || and the earlier point-coupling results of Ref. [15|. The 
larger compressibilities of the point-coupling models is evident and leads to deviations from 
the meson-model result. With increasing density, the energies for point-coupling models 
VA4 and FZ4 (and all others fit here) increase smoothly, as found for meson models (such 



as Gl). In contrast, the point-coupling model found by Nikolaus et al. |18fl shows a peculiar 



second minimum. This structure comes from the contributions from the highest-order terms 



(fourth power of density) in the model of Ref. [|ll| , which are much larger than those of the 
next-to- leading order. Their contributions to the value of M* /M and to the binding energy 
become dominant for values of the Fermi momentum not far from saturation, indicating 
that the expansion has broken down; for this restricted class of models the extrapolation in 
density is not valid. The corresponding parameters are individually unnatural and there are 
no mixed terms, so the M* equation is strongly affected. 

The implications of naturalness are clear in Figs. |] and [|, where contributions to the 
binding energy per nucleon at equilibrium from terms of the form p % s pP v with % + j = n and 
p v = j® are plotted for various parameter sets. (The reader is cautioned not to confuse 
n with the label of the parameter sets associated with the optimal hierarchy). The scale 
of the natural size of contributions expected from an n-th order term is obtained through 
NDA; since the difference in the scalar and vector densities is small, (the magnitude of) the 
contribution from such a term can be estimated by 

cwv(j^f\ (52) 
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FIG. 4. Contributions to the binding energy of equilibrium nuclear matter from terms of the 
form p\pi u with i + j = n. The symbols indicate terms with i = (squares), j = (circles), and 
^ 7^ 0, j (diamonds). The unfilled symbols indicate set FZ4 and the filled symbols indicate set 
FA4. The crosses are estimates based on Eq. (|5^), 
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FIG. 5. Contributions to the binding energy of equilibrium nuclear matter from terms of the 
form plpl with i + j = n. The symbols indicate terms with i = (squares), j = (circles), and 
i 7^ 0, j / (diamonds). The unfilled symbols indicate set VZ4 and the filled symbols indicate set 
VA4. The crosses are estimates based on Eq. (^). 
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where a nat is a (positive) number of order unity and p v is the nuclear matter equilibrium 
density. Estimates for A = 770 MeV are marked in the figures with crosses and the error 
bars reflect a range of a nat from 1/2 to 2. We observe that the estimates are quite consistent 
with the energies found from the fits. 

Since the contributions decrease steadily with increasing powers of the densities, the 
expansion and truncation scheme proposed under NDA is justified. Note that this decrease 
will become more gradual as the density is increased above equilibrium density. Based on 
the estimates, a truncation at n = 4 for ordinary nuclear observables yields an error of order 
1 MeV or less. The contributions at fourth-order in the density are on the order of the 
nuclear-matter binding energy, however, and consistency implies we need to include these 
terms. In practice, however, the n = 3 truncation (parameter set P2) is able to absorb the 
n = 4 contributions. This result is consistent with a recent analysis of the nonrelativistic 



Skyrme energy functional ||31|| . On the other hand, this was not the case in Ref. ||, where 
the n = 3 fit, while quite reasonable, was also noticeably inferior to n = 4. 

While the expansion and truncation scheme is supported by our results, a comparison of 
the various parameter sets shows significant variation for individual parameters. In many 
cases even the sign is indeterminate. Thus the parameters in this approach, as in the 
relativistic meson model, are underdetermined. On the other hand, Table [VI] shows that the 
optimal parameters are much better constrained by the data. Therefore reformulating the 
effective theory in terms of an optimal basis may be more productive. 



C. Equivalent meson masses and Yukawa couplings 

No restrictions were imposed during optimization on the signs of the coefficients k,2, 
(2, and £2 nor on Kd, Q, and £<j when they were allowed to vary (the "V" sets). On the 
other hand, as shown in Tables [IlHV|, a direct transformation from a relativistic mean-field 
meson model would predict definite signs for these coefficients; the coefficients determined 
through this transformation would depend only on the squares of the Yukawa couplings, 
gf, as well as the squares of the meson masses, mf. Sets FAn and FZn all have positive 
values for the corresponding values of gf as well as for the scalar mass squared. (Recall 
that the corresponding values of m v and m p were held fixed for these sets). Parameters sets 
VZn and VAn, however, include sets with mf < and gf < 0. The ratio gf/m?, which is 
the integrated strength in each channel, is reasonably well determined and is in every case 
positive. 



The difficulty with predicting individual masses and couplings is revealed by Table [VI . 
While Ai is well determined and consistent across the parameter sets, A2 is poorly deter- 
mined. Thus, the individual values for R d and Q are also poorly determined, and we cannot 
reliably extract masses and couplings from the fits. We conclude that the nuclear observ- 
ables we have used do not provide sufficient constraints to definitively test whether the point 
couplings are dominated by an underlying meson phenomenology. The combinations that 
are well determined correspond to the strength and range of the effective central potential. 
Specifically, if the static nonrelativistic potential in momentum space from the exchange of 
scalar and vector mesons is written as: 

^(q) = 47r (Sr-5 + -&-^) ' ( 53 ) 

V mr s + + q / 
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and expanded in powers of q 



v(q) = ^ f-4 + 41 - 4 ^q 2 f-4 + 41 + • • • , (54) 



ml m„ 2 I \ mi mi 



the two combinations of couplings and masses in parentheses are reasonably well determined. 



VII. SUMMARY 

The successful application of EFT concepts to mean-field mesonic models of nuclei || 
motivates their application to point-coupling models. Such models describe nuclear inter- 
actions through contact terms (in a derivative expansion) in place of meson exchange. In 
comparison to the large body of work on meson models, the analysis of relativistic point- 
coupling models has been quite limited. Here, we extend the model of Nikolaus, Hoch, and 



Madland |L8[ and its analysis by Friar et al. |L9| to encompass a more complete analysis 
based on EFT concepts. 

The lagrangian is consistent with the symmetries of QCD and is organized according 
to the same principle of naive dimensional analysis (NDA) applied to the relativistic meson 
model of Ref. ||. The organization is an expansion in powers and derivatives of the densities 
in ratio to scales dictated by NDA results: a given term in the lagrangian takes the form 

p (NTN\ 1 ' 



(55) 



where c is a dimensionless coefficient and p and I are integers. (Electromagnetic interactions 
will also contain a power of the photon field and the electric charge). In principle, all terms 
consistent with the symmetries should be included to a given order in the expansion, but in 
the present work we have omitted a variety of terms that we expect to be poorly determined 
based on meson exchange phenomenology and experience with relativistic meson mean-field 
models. 

The NDA organization provides a valid expansion scheme provided the coefficients are 
"natural" (on the order of unity); one can construct an energy functional in powers of 
densities and their derivatives and truncate at some finite order with a controlled error. Our 
fits to bulk nuclear properties show this to be the case. Beyond a second-order truncation, 
the models resulting from these fits reproduce the experimental data quite well. We conclude 
that NDA and the naturalness assumption are compatible with and implied by the observed 
properties of finite nuclei, even though many-body effects are absorbed into the coefficients. 

Point-coupling models therefore provide an alternative phenomenology to mean-field me- 
son models. In principle, a direct transformation exists between any point-coupling and 
meson mean-field model. The equations of motion for the mean meson fields can be iterated 
to solve the fields as an expansion in scalar and vector densities. The delicate cancellations 
at low order in the expansion result in too large a truncation error for such transformations 
to be of use, however, and any point- coupling model derived from a mesonic model in this 
way must be refit to the data. We note that naturalness of the point-coupling model was 
found in the absence of counting factors that were required in the meson model analysis. 
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Although fourth-order terms would be required for a consistent truncation, a truncation 
at third order in the point-coupling model can yield a good fit to data, in contrast to the 
meson models where fourth-order terms were necessary in obtaining a good fit. Due to im- 
provements in the optimization procedure over that used in Ref. @ through a reorganization 
of the lagrangian in terms of "optimal parameters" , a true comparison cannot be made until 
the fits obtained in that reference are reanalyzed. 

The analysis in terms of optimal parameters suggests that they may provide a more effi- 
cient basis for an expansion. While the individual coefficients in the covariant point-coupling 
model are in general poorly determined by the nuclear data, the lower-order optimal pa- 
rameters are quite well constrained and only the highest-order parameters are badly under- 
determined. The use of optimal parameters was suggested by methods of the heavy baryon 
formulation of chiral perturbation theory, adapted to finite density. Further development of 
this approach and its relation to nonrelativistic energy functionals for nuclei is in progress. 
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